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Abstract 

We show that the best fits of the MSSM constrained by a simple SO(10) SUSY 
GUT are consistent with the present data on the muon anomalous magnetic moment. 
The best fits assume rather large values of tan/3 ~ 50, and in our analysis they are not 
a priori correlated in any way, directly or indirectly, with the experimental limit on 
a M . Regions in the SUSY parameter space, which are currently ruled out because of 
too large a^ USY , are already excluded in the global fit by excessive corrections to mj, 
unacceptable BR(b — ► S7), or direct experimental limits on sparticle masses. However, 
our results indicate that the accuracy expected in the ongoing E821 experiment at BNL 
will eventually turn the muon anomalous magnetic moment into a major constraint for 
this regime of the MSSM. 
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1 Introduction 



The anomalous magnetic moment of the muon, = (g^ — 2)/2, is potentially a significant 
constraint for any extension of the Standard Model (SM). The measured value quoted by 
the Particle Data Group is 

a^ DG = (11 659 230 ±84) x KT 10 . (1) 

This value does not take into account the early results [] from the E821 experiment at BNL 



a BNL97 = (11 659 250 ±150) x KT 10 , (2) 
a BNL9s = (H659191 ± 59) x 10 -io. (3) 

On the other hand, the SM prediction yields (0 and references therein) 

af* = (11 659 159.6 ±6.7) x 10" 10 , (4) 

where the QED, electroweak and hadronic contributions are sumed and the errors are com- 
bined in quadrature. When we take into account all these results, contributions from new 
physics beyond the SM are constrained to fit within the window 

- 28 x 10- 10 < 5a^ EW < ±124 x 10' 10 , at 90% C.L. (5) 

The width of the window is at present dominated by the experimental uncertainty. That, 
however, will change after the E821 experiment is completed. The inclusion of the current 
BNL data already reduces the available window by a factor of 2. After two more years of 
running the E821 is eventually expected to reach the accuracy ±3 x 10~ 10 . If the measured 
central value then turns out to be exactly equal to the SM prediction, eq.(£|), the new 
constraint will be 

- 12 x 10" 10 < 5a^ EW < +12 x 10" 10 , at 90% C.L., (6) 

with the window for new physics narrowed by more than a factor of 6. 

In the MSSM, there are significant contributions from new physics due to the chargino- 
sneutrino and neutralino-smuon loops: 
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1 The 1998 result is still preliminary. 

2 This estimate does not take into account possible improvements over time in the hadronic uncertainty, 
and thus the constraint on new physics may actually be even tighter. 
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and the standard integrals over the two-dimensional Feynman parameter space are 



In these relations, notation of ||, ||, |7j is assumed. C"s and iV's are the couplings of the 
chiral muon states with the charginos and neutralinos, respectively. Uja, Vai and N a3 are 
the elements of the chargino and neutralino mixing matrices, A = 1, 2; a = 1, 4. T's are the 
mixing matrices of the sleptons after the flavor eigenstates have been first rotated by the 
unitary matrices which diagonalize the respective fermionic states. Note that a = 1, 3 for 
sneutrinos and a = 1, 6 for charged sleptons. In particular, for the charged sleptons each 
(r e ) ai is a 6x3 matrix defined as (T e L) ai = Z ai , (T e ji) a i = Z ai+3 , (i = 1, 3), where Z is a 6x6 
mixing matrix for charged sleptons of all three generations. For 3 sneutrinos, T v l directly 
diagonalizes the 3x3 sneutrino mass matrix. m^+, m v o, m~ and m~ are the chargino, 
neutralino, sneutrino and charged slepton mass eigenvalues, m M and A M are the muon mass 
and diagonal muon Yukawa coupling, and g\ and #2 are the electroweak gauge couplings. 

Note that at first glance the terms proportional to F 3 and F 4 in eqs.(^, §), respectively, 
are enhanced by m x /m^. The net enahncement is actually of the order v u /vd = tan/5, (va 
and v u are the Higgs vevs which give masses to the <i-quarks and charged leptons, and to 
the li-quarks, respectively.) [] The enhancement by tan/3 can be traced back to the diagrams 
with the chirality flip inside the loop (or in one of the vertices) as opposed to the terms 
proportional to F% and F 2 where the chirality flip takes place in the external muon leg. 
There are no similarly enhanced terms in the SM, where the chirality can only be flipped 

3 Note that there is no additional enhancement due to 1/ A M in the terms proportional to F3 and F4 after 
the sum over the mass eigenstates is performed. In fact, any combination of chiral muon states contributes 
to 5a,p with the net contribution suppressed by small Yukawa coupling A M , as one can see, for instance, 
from the relevant Feynman diagrams in terms of flavor eigenstates (i.e., in the interaction basis). That 
also explains why the electron anomalous magnetic moment is less sensitive to new physics than its muon 
analogy. 




(10) 
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in the external muon leg. Thus for tan/3 > 1 we expect that the terms proportional to F3 
and F 4 dominate in eqs.(|7|) and and that compared to the SM electroweak gauge boson 
contribution Sa^ w , the SUSY contribution scales approximately as 

6a s u u SY =< 5a E u w ( ^] tan/3 ~ 15 x 1(T 10 (l°22fV\ 2 tan/3; (n) 

where m stands for the heaviest sparticle mass in the loop. The effect has been noticed and 
emphasized by the earlier studies focusing on the muon anomalous magnetic moment in the 
context of the MSSM |§. 

Relation ( pT|) predicts a very simple tan/3 dependence. It suggests that the large tan/3 
regime of the MSSM may already be constrained by the currently allowed window, eq.(|5|). 
That is of interest for models of grand unification based on simple SO(10), since the SO(10) 
GUT constraint yt = Vb = Vt f° r the Yukawa couplings of the third generation implies 
tan/? 50. In this study, we first apply our analysis to the MSSM constrained only by gauge 
coupling unification and, as a warm-up, compute the muon anomalous magnetic moment for 
fixed values tan/3 = 2 and tan/3 = 20. Next, we proceed to the best fits of a simple SO(10) 
SUSY GUT ||. We present our results in terms of the contour lines of constant u drawn 
in the SUSY parameter space. In section 2, we review our numerical analysis. Sections 3 
and 4 discuss the results and prospects for the future. 

2 Numerical Analysis 

Our numerical analysis has relied on the top down global analysis introduced in |J. The soft 
SUSY breaking mediated by supergravity (SUGRA)|TI| was assumed throughout this paper. 
As explained at the end of the previous section, the magnitude of tan/? is of special interest 
for the analysis. Hence we discuss separately three cases with low, medium, and large value 
of tan/?. First, fixed tan/3 = 2 was considered, along with the scalar trilinear parameter A 
set to zero to simplify the analysis which is insensitive to Aq in this case. Next, tan/? was 
raised to 20 and A set free to vary. In the third case, we used the results of the global 
analysis of model 4c, a simple SO(10) model with minimal number of effective operators 
leading to realistic Yukawa matrices. The model was suggested by Lucas and Raby || and 
its low-energy analysis was presented in ||. A direct model dependence of a M on the details 
of the Yukawa matrices is, however, very limited. In this case, tan/? was not fixed to any 
particular value. Instead, it was a free parameter of the global analysis but since the model 
predicted exact t — b — r Yukawa coupling unification, the best fit values of tan/? were always 
found between 50 and 55, dependent on a particular (jno,Mi/2) point as explained below. 
Aq was free to vary — as for tan/3 = 20 — which allowed the best fits to optimize the effects 
of the left-right stop mixing, which is enhanced by tan/?, in the analysis. 

In each of the three cases, gauge coupling unification was imposed up to a small (less 
than 5%) negative correction, called 63, to a s at the unification scale Ma- Scale Mq has 
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been denned as the scale where ot\ and «2 are exactly equal to the common value ocg- For 
low and medium tan/3 (2 and 20), the minimal set of the initial SUSY parameters 



M1/2, mo, A , sign/x, and tan/3 (12) 

was assumed, with M 1( / 2 , mo, and Aq (the universal gaugino mass, scalar mass and trilinear 
coupling) introduced at Mq. The Yukawa couplings of the third generation fermions at Mq 
were unconstrained and free to vary independently on each other. For large tan/3 in the third 
case, the scalar Higgs masses were allowed to deviate from mo in order to alleviate strong 
fine-tuning required for the correct electroweak symmetry breaking. Thus instead of fll2|) , 
the set 

Mi /2 , mo, (m H Jm ) 2 , (m Hu /m ) 2 , A , sign/z, and tan/3 (13) 

was actually used as initial SUSY parameters. As already mentioned, the third generation 
yukawas were strictly set equal to each other at M G in this case. In fact, the SO(10)-based 
equality among them is the main reason why such a large tan/3 is attractive. 

The rest of the analysis was then the same for each of the three cases. Particular values 
of mo and M1/2 were picked up, while the rest of the initial parameters varied. That included 
varying ac, Mq, 63, the third generation Yukawa couplings, and A (for medium and large 
tan/3). Using the 2-loop RGEs for the dimensionless couplings and 1-loop RGEs for the 
dimensionful couplings the theory was renormalized down to the SUSY scale, which was 
set equal to the mass of the Z boson. The electroweak symmetry breaking was checked 
to one loop as in Q, based on the effective potential method of ref . [jrTJ . One- loop SUSY 



threshold corrections to fermion masses were calculated consistently at this scale. That is 
of particular importance for mj which receives significant corrections if tan/3 is large. Also, 
the experimental constraints imposed by the observed branching ratio BR(b 37) and 
by direct sparticle searches were taken into account. Finally, 5a^ USY = 5a* + + 5a%° was 
evaluated following eqs. (0) and @ for those values of the initial parameters which gave the 
lowest x 2 calculated out of the ten low energy observables M z , M w , p new , a s (M z ), a, G^, 
M t , m b (M b ), M T , and BR{b — > 57). Details of the low energy analysis can be found in 
0. The calculated value of 5a^ USY did not have any effect on the \ 2 calculation and the 
subsequent selection of the SUSY parameter space. 



3 Results and Discussion 

The results for tan/3 = 2 and tan/3 = 20 are shown in figures [l|-44|. The figures p] and |2] show 
the contour lines of constant 5a^ USY x 10 10 in the (mo, M1/2) plane. These results are then 
transformed into the dependence on physical masses in figures || and |], where we chose the 
(m£>,m x +) plane, with m^ being the muon sneutrino mass and m x + being the mass of the 
lighter chargino. The contour lines in these figures are bound from below by the limit on the 
neutralino mass (a limit m x o > 55GeV was imposed flatly, for simplicity), and from above 
by the stau mass (m? > 60GeV, in this analysis). The important observation is that the 



4 



present limits on Sa^ USY , eq.(g|) are not excluding any region of the parameter space which 
is left available by other experiments. The size of <5a^ c/sy is in agreement with the results of 



Chattopadhyay and Nath, and of Moroi ||. We can also confirm their observation that the 
neutralino contribution <5a*° is generally much smaller than the chargino contribution 5a* + in 
the whole parameter space available. Yet several interesting characteristics of the presented 
results cannot be read off directly from the earlier studies. It is first of all the simple pattern 
suggested by figures | and £|. The pattern in the figures shows a strong dependence of 6a^ USY 
on the sneutrino mass and very little sensitivity to the mass of the lighter chargino. That 
can be understood from the fact that the loop integrals are in general most sensitive to the 
heaviest mass in the loop, and the universal initial conditions (see (0)) together with the 
experimental limits from direct searches always lead to > m x + at the electroweak scale. 
Actually, the figures show how well the approximate relation (|ll|) works. One could e.g., 
directly read out the value of sneutrino mass from a figure of this type once a more accurate 
measurement of a M and the value of tan/? become available. Alternatively, a more accurate 
measurement of a M can be converted into a limit on tan/3 in a straightforward way, if the 
sneutrino mass is known. Also note how well the linear dependence 5a^ USY cx tan/3 holds. 
If we overlapped figures |3| and f|, the contour lines marked as 10, 5, 2, 1, and 0.5 in figure |3] 
would be practically on top of the contour lines marked as 100, 50, 20, 10, and 5 in figure [|. 
That suggests that the tan/3 enhanced terms in eqs. (0) and (|8|) become dominant already for 
tan/3 > 2. We checked this feature in the numerical analysis: typically, the term proportional 
to F 3 in eq.([7D accounts for about 80-85% of the net 5a* + for tan/? = 2, while it is 97-99% 
for tan/3 = 20. (The analogous term in eq.@ is less dominant, but has a smaller net effect 
since Sa}° < Sa}^ .) 

From these results one could extrapolate that the current a M limit, eq.([5]), places an 
important constraint on the analysis if tan/3 is as large as 50. 

However, case tan/3 ~ 50 is qualitatively different. As discussed in the study by Blazek 
and Raby [Q, the global analysis yields two distinct fits, figs. |5|a b. The fits differ by 
the sign of the Wilson coefficient Cj ISSM in the effective quark Hamiltonian below the 
electroweak scale. (The coefficient CV determines the b — > 37 decay amplitude, after the 
QCD renormalization effects are taken into account (T^j.) In the MSSM with large tan/3, the 
sign of C7 can be the same, or the opposite, as in the SM. It can be reversed due to the fact 
that the chargino contribution can be of either sign. This contribution is enhanced by tan/? 
compared to the SM and charged Higgs contributions, whose signs are fixed and alike, and 
thus the flipped sign of Cf 1831 ^ 1 cannot be obtained for low tan/3. [] For tan/3 ~ 50 the fit 
with the flipped sign is equally good as the fit with the sign unchanged, see figure ^. The fits 
differ in the range of the allowed SUSY parameter space: to reverse the sign the chargino 
contribution accepts lower squark masses and different values of A than in the case when 
the sign is unchanged. 

One can anticipate that these differences will be reflected in the analysis of <5a^ c/5y with 



The contributions with ncutralinos or gluinos in the loop always turn out to be small in our analysis. 
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tan/3 ~ 50, since the quark sector is correlated with the lepton sector through the unification 
constraint, and the soft masses are interrelated through the universality assumption. Due 
to differences in SUSY spectra we can expect Sa^ USY to be more significant in the fit with 
the flipped sign of C^ SSM than in the fit where the signs are the same. These expectations 
are indeed realized in our results in figs. |a and |6|b. In these figures, we plot the contour 
lines in the (m , Mi/ 2 ) plane of 5a^ USY , calculated according to eqs. (0) and (§), with 
all the masses, mixings and couplings taken over from the best fit values at the specific 
(mo, M1/2) point in the SUSY space.f] The important feature, which runs contrary to the 
naive extrapolation from the previous two cases, is that the MSSM contribution to the muon 
anomalous magnetic moment stays within the currently allowed range at 90% C.L., even for 
tan(3 as large as predicted by simple SO(10) GUTs. Despite the fact that figures |l]-|4] suggest 
that the analysis with tan/3 = 50 — 55 becomes sensitive to the current constraint on a M , 
eq.(|5]), our results in fig.[| clearly indicate that this is not the case. The reson for this is that 
we consistently demand that all known particle physics constraints (besides those imposed by 
a^) are accounted for. For tan/3 as large as 50, the allowed SUSY parameter range is further 
reduced by strong constraints on the b quark mass and the branching ratio BR(b — > 57), 
when compared to the regime where tan/3 = 2 or 20. 

Finally, we make a note on the sign of fi. As clearly indicated on the top of each 
figure, we have presented our results just for /i > 0. For this sign of the /i parameter, the 
SUSY contributions to a M are positive across the whole (mo,M 1 / 2 ) plane. For /i < 0, the 
contributions change sign too. However, the chargino contribution to C7 also changes sign 
in this case, which leads to unacceptable values of BR(b — > 57) for medium and large tan/3. 
It is amazing to observe that for this range of tan/3 the sign of the [i parameter favored by 
b — > 57 is the same as the sign preferred by the experimental window open for 5a^ EW . R 



4 Prospects for the New BNL Experiment and Con- 
clusions 

When the future BNL experiment reduces the uncertainty on a M down to ±12 x 10~ 10 
at 90%C.L., eq.(^|), the muon anomalous magnetic moment will undoubtfully turn into a 
powerful constraint on the MSSM analysis. It is already clear from figures |I|, |2|, and |5] that it 
will be a major constraint for large and medium tan/3 in the region mo < 400 — 500GeV. Of 

5 We do not show plots in the (mc,,m x +) plane in this case. For large tan/3 the best fits of the global 
analysis result in the lightest chargino being higgsino-like, with its mass close to the electroweak scale across 
the whole (moj-M^) plane. Thus the two-dimensional plots in the (top,to x +) plane would contract to a 
dependence on mo only. 

6 For low tan/3, the b — > 57 constraint disappears (the chargino contribution to C7 is small) and \i can 
be negative. However, figure |l| shows that 8a^ USY is very small in this case. We know from the earlier 
studies Q that no substantial change in the magnitude of Sa^ USY occurs for different signs of Thus we 
conclude that for fi < only the low tan/3 case is viable, but then the SUSY contributions to a M stay below 
the electroweak contributions and will be hard to observe even after the E821 experiment is completed. 
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course, the BNL result may drastically affect the MSSM analysis with any value of tan/3 if the 
observed central value turns out to be below (or well above) the SM prediction. For such an 
outcome, the muon anomalous magnetic moment will actually become a dominant constraint 
for the MSSM analysis with the universal SUGRA-mediated SUSY breaking terms. 

In the meantime, a M does not pose any new constraints on the MSSM analysis under the 
terms considered in this work. 
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Figure 1: Contour lines of constant 5a^ x 10 in the analysis with fixed tan/? = 2 in 
the SUSY parameter space. 
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Figure 2: Contour lines of constant Sa^ USY x 10 10 in the analysis with fixed tan/3 = 20, in 
the SUSY parameter space. 
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Figure 3: Contour lines of constant 5a^ x 10 10 in the analysis with fixed tan/5 = 2, 
in the (m f ,m x +) plane, where is the muon sneutrino mass and m x + is the mass of the 
lighter chargino. 
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Figure 4: Contour lines of constant 5a^ x 10 10 in the analysis with fixed tan/3 = 20, 
in the (m f ,m x +) plane, where is the muon sneutrino mass and m x + is the mass of the 
lighter chargino. 
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Figure 5: x 2 contour plots in the best fits of a simple SO(10) model, with the Wilson 
coefficient C^ SSM (relevant for the b — > S7 decay) of (a) the same (b) the opposite sign as 
compared to Cf M . As indicated, the contour lines correspond to x 2 — 6, 3, and 0.3 per 
3 d.o.f., respectively, tan/3 varies between 50 and 55 due to the model prediction that t, b, 
and r Yukawa couplings are equal at the GUT scale. 
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Figure 6: Contour lines of constant 5a^ USY x 10 10 in the best fits of a simple SO(10) model, 
with the Wilson coefficient C^ SSM of (a) the same (b) the opposite sign as compared to . 
For better reference, the x 2 contour lines of figures 5a and 5b are shown in the background 
of (a) and (b), respectively, as dotted lines. 
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